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Abstract. - A Landau theory is constructed for the vulcanization transition in cross-linked 
polymer systems with spontaneous nematic ordering. The neo-classical theory of the elasticity 
of nematic elastomers is derived via the minimization of this Landau free energy; this neo- 
classical theory contains the classical theory of rubber elasticity as its isotropic limit. Our 
work not only reveals the statistical-mechanical roots of these elasticity theories, but also 
demonstrates that they are applicable to a wide class of random solids. It also constitutes a 
starting-point for the investigation of sample-to-sample fluctuations in various forms of vulcan- 
ized matter. 



The classical theory of rubber elasticity [1] has been remarkably successful. A blend of 
phenomenology and molecular-level reasoning, it is based on a few simple assumptions, and 
bears great predictive and descriptive power. By modeling rubbery materials (i.e. elastomers) 
as incompressible networks of entropic Gaussian chains, it gives their elastic free energy density 
/ at temperature T as 

/ = ^TrA T A, (1) 

for all uniform deformations r — > A • r that conserve the volume (i.e. det A = 1). For 
most rubber-like materials the assumption of incompressibility is well satisfied^ 1 ) The shear 
modulus fi is given by n c T, where the constant n c is usually referred as "the density of effective 
chains in the network'.' The classical theory [i.e. Eq. (JTJ and associated arguments] explain 
many essential features of rubbery materials, such as their stress-strain curves (at least for 
deformations that are not too large), and the striking temperature dependence of their shear 
moduli, as well as their strain-birefringence (i.e. the stress- optical effect). 



( 1 )This assumption may break down for swollen rubbers and gels. In this case, a finite-compressibility version 
of Eq. {TJ remains valid. In fact, this is what we will derive in the paper. 
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There are several important issues unresolved by the classical theory. First, for a given 
cross-link density the shear modulus is not calculated within the theory. ( 2 ) Second, in the 
intermediate-strain range there is universal and significant downward deviation of the exper- 
imental stress-strain curve, compared with the theoretical prediction. Finally, the issue of 
polymer entanglement is not addressed by the classical theory. 

Subsequent efforts to improve the classical theory of rubber elasticity have focused on vari- 
ous directions (for one overview see Ref. [2]). The non-Gaussian nature of the chain statistics, 
due to the finite extensibility of the polymers, has been taken into account, and explains the 
large-deformation behavior of the stress-strain curve. Purely mathematical modeling, as in 
the theories due to Mooney, Rivlin and others (see, e.g., Ref. [1]), also provided useful insight. 
At the microscopic level, the effects of chain entanglement have long been emphasized and 
modeled via various approximation schemes, most notably by the Edwards tube model [4] 
and its derivatives, although results from these models are often either inconclusive or contra- 
dictory. It seems fair to say that none of these efforts is as successful as the classical theory, 
either in terms of simplicity of assumptions, or in terms of broad descriptive power. 

Recently, a simple anisotropic generalization of the classical model [5,6], known as the 
neo-classical model, was constructed to describe the highly unusual elasticity of nematic elas- 
tomers, i.e., rubbery materials with (spontaneously) broken rotational symmetry, and has 
done so with great success. In the presence of nematic order, the "step-length tensor" 1 char- 
acterizing the conformations of the polymer chains is anisotropic. Then, according to the 
neo-classical model, the elastic free energy of a nematic elastomer with deformation A is given 
by 

/ = |Trl A T r 1 A, (2) 

where 1q and 1 are the (in general anisotropic) step-length tensors in the initial and deformed 
states, and are functions of the nematic order parameter tensors Q° and Q in these two 
states, respectively. A remarkable feature of Eq. is that for a given 1q and 1 there exists 
a continuous manifold of deformations A that cost zero elastic free energy. These novel soft 
modes [7, 8] , as well as the thermal and quenched fluctuations associated with them, have been 
the focus of intensive study, both experimentally (see, e.g., Ref. [6]) and theoretically [6,9,10]. 

In a classic paper, Deam and Edward [11] initiated a statistical- mechanical approach to 
the study of rubber elasticity that incorporates both thermal and quenched fluctuations, along 
with repulsive interactions. This replica-based approach, which has been called "vulcanization 
theory" (VT) [12, 13], has been explored in great detail. Progresses in this direction includes, 
inter alia: (i) the calculation of the mean-field order parameter [12]; (ii) the derivation of a 
universal Landau theory [14]; (iii) the development of connections with percolation theory [15]; 
and (iv) the analysis of critical scaling for shear modulus [16]. A main virtue of VT is that it 
follows the Landau paradigm of modern condensed matter physics, inasmuch as it concentrates 
on order parameters, symmetries and length-scales. In particular, because of its independence 
on microscopic details, we view the Landau theory of the vulcanization transition as the right 
theory if one wishes to address the long length-scale physics of rubbery materials, especially 
near the vulcanization transition. 

The aim of the present work is to establish connections between VT and the classical — as 
well as the neo-classical — elasticity of isotropic and nematic elastomers. We begin this task by 
generalizing the Landau theory for the vulcanization transition to systems with spontaneous 

( 2 )We remind the reader that the classical theory was developed before percolation theory was. Therefore, 
how an infinite network emerges during the random cross-linking process was not understood. In fact, near 
the vulcanization point, the "effective chains" of the classical theory bear little resemblance to the original 
polymer chains before cross-linking. 
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nematic ordering. As a saddle-point approximation to this theory, we derive the neo-classical 
model for the elasticity of nematic elastomers. In the isotropic limit, this result reduces 
to the classical theory of rubber elasticity Q). Our work not only reveals the statistical- 
mechanical roots of these elasticity theories, but also demonstrates that they are applicable 
to a wide class of random solids. It also constitutes a starting-point for the investigation of 
sample-to-sample fluctuations in various forms of vulcanized matter. 

We begin with the real-space version of the order parameter for the replica field theory of 
the VT [12,14,16,17] in d dimensions, which is a function of (1+n) (i-vectors x = (x°, . . . ,x n ): 

N n 

n(*) = n(*f\ -..,*») = £ ( n - c " )> 1+ „ - ™ ■ (3) 

j=l a=0 U 

Here, c" (with a = 0, 1, . . . ,n) are the 1 + n replicas of the position d-vectors of the N 
monomers (with j — 1, . . . , N) that comprise the system. Vq is the volume of the system in 
the preparation state, and V is the volume of the system in the measurement state (which 
may differ from Vo). The brackets (■ • • }i + „ denotes an average over the replica field theory 
given below. Up to a constant, the order parameter of VT gives the conditional probability 
that a monomer found at x° at the time of cross-linking is later found at {x 1 , . . . ,x™} in n 
independent measurements after cross-linking. The one-replica parts of f2 (for a = 0, . . . , n) 
are defined via 

0>(x Q )=/ J] dx*n(S) = X;(5(x«-c?)) (4) 
J i=i 1+n Va 

where V a = (Vo, V, .. . , V). Of these, S!o(x°) describes the density fluctuations in the prepara- 
tion ensemble, and 17 Q (x a ) (for a = 1, . . . ,n) describe density fluctuations in the measurement 
ensemble. 

The Landau free energy for VT for isotropic systems is given by 

H V [Q} = |di j^( Vo fi) 2 + f g(v Q n) 2 + In 2 -|n 3 | (5) 
+f /dxW + f ±Jd^(n a f + ^ +Po v + ^ + n P v. 

a — 1 

(pO) Bo, Ko) an d (p, B, K) are, respectively, the pressure, inverse susceptibility for density fluc- 
tuations and chain stretchability in the preparation and measurement ensembles. The chain 
stretchability is proportional to the squared radius of gyration of each constituent polymer 
chain in the isotropic state. In the absence of any externally imposed deformation, the values 
of Vo and V should be such that they minimize the Landau free energy. For simplicity, we 
consider the case of equal pressures and bulk moduli in the preparation and measurement 
ensembles (i.e. po — p and Bq = B). Near the vulcanization transition, where the order 
parameter ft is small, this leads to V 2 — V 2 — B N 2 /2p, in the absence of deformation. 
The control parameter r triggers the transitions to the solid phase when it becomes nega- 
tive (i.e. when the density of cross-links exceeds some critical value). This model has been 
analyzed extensively, both within and beyond the mean-field level( 3 ). Last but not least, 

( 3 )In the early literature, the model was usually expressed in momentum space rather than of real space. 
Furthermore, the one-replica sector of Q was excluded explicitly, which amounts to setting Bo and B to infinity 
in Eq. JSJ- The asymmetry between the preparation and measurement ensembles, and its consequences, have 
been stressed only recently [16,18]. 
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we note that the order parameter 12 (x) of this Landau theory is a single-monomer quantity 
(albeit replicated). The original polymer degrees of freedom are integrated out in deriving the 
Landau theory and, consequently the issue of topological entanglement becomes irrelevant in 
this theory. It is our belief that inclusion of entanglement in the original theory would simply 
lead to a quantitative modification of the parameters in the Landau theory Eq. (JjjJ, not an 
invalidation of the theory itself. 

There are various ways to incorporate nematic ordering into VT. The simplest is to couple 
the VT order parameter 12 to the replicated symmetric traceless tensor fields Q Q (with a = 
0,1, ... ,n). Of these, Q° describes the nematic order in the preparation ensemble, whereas Q" 
(for a = 1, ■ ■ • ,n) describe nematic order in the (n-fold replicated) measurement ensemble. 
The resulting free energy must be invariant under the simultaneous rotation of Q"and the 
spatial position vectors x a , independently for each replica a. The lowest-order coupling (in 
Q and gradients) allowed by symmetry is 

Hn Q = Jdx f^QSbVSnvgn + ^f^Q^vjnvfn^ , (6) 

where V" indicates a derivative with respect to the a th cartesian component of the a th repli- 
cated position vector. The signs of the coupling constants 770 and 77 depend on details of the 
chemical structure of nematic polymers under consideration. In a separate publication [19] 
we shall derive Eq. © from a lower-level description 

The total free energy should also be augmented by a part that depends only on Q Q , 
and accounts for the nematic interactions between neighboring (anisotropic) monomers. It 
is known, however, that the nematic energy-scale (roughly k^T per monomer) is orders of 
magnitude larger than the energy-scale for the vulcanization transition (roughly k-gT per 
chain). Therefore, we may neglect the feedback of the 12 ordering on the nematic order, 
and thus may treat Q a as given. More specifically at the mean-field level, and under the 
assumption of that there is neither replica nor macroscopic translational symmetry breaking, 
we may set both Q° and Q a (= Q for a = 1, . . . , n) to be constants, characterizing the 
uniform nematic order in both the preparation and measurement states. Additionally, large 
and positive values for Bo and B guarantee that the saddle-point value of 12 vanishes in the 
one-replica sector, i.e. Eq. vanishes at the saddle point. By minimizing the total free 
energy over the VT order parameter 12, we find the saddle-point equation 

= -i° ab v° v° n - UK v?^ + r n - i«n 2 . (7) 

a=l 

As they stand, the tensors lo and 1 are short-hand for 

l° ab = K S ab + 77q Q° ab , l ab = K5 ab + riQ ab . (8) 

However, as we shall see below, they are in fact the effective step-length tensors in the prepa- 
ration and measurement ensembles that appear in the neo-classical elastic free energy, Eq. • 
In the absence of any externally imposed deformation, the saddle-point equation J7J) is 
solved by the following Ansatz( 4 ): 




V 



l+n 



,(9a) 



N{t) = (tt/t) 11 ^ (detio)* (detl)^, (9b) 



( 4 )A very similar Ansatz in real space has been used in Reference [18] for the isotropic case. 
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a=0 



a=l 



OC=ra : 



Fig. 1 - An externally imposed deformation changes afnnely the boundary of the system in the n 
replicas of the measurement ensemble, but not that of the preparation ensemble, i.e. the th replica. 



where the <i-dimensional vector z is integrated over the interior of system Vq in the prepara- 
tion state. Obviously, if both the preparation state and the measurement state are isotropic 
(i.e. Q° = Q = and lo oc 1 oc I, the above saddle-point Ansatz reduces to the earlier form 
appropriate to isotropic systems [12], with p(r) the so-called distribution of inverse square 
localization lengths. The interpretation of this saddle point is as follows. A certain fraction 
of the monomers (q per unit volume) belong to the infinite cluster (i.e. the gel fraction) and 
are localized. Such a monomer fluctuates around the point z with Gaussian variance-matrix 
t _1 1 in the preparation ensemble (i.e. the th replica) and fluctuates around the same point 
z( 5 ) with variance-matrix r _1 l in the measurement ensemble (i.e. replicas 1 to n). From 
Eq. JSJ it is easy to see that the role of non- vanishing nematic order is to confer anisotropy on 
these fluctuations. Finally, because of the random nature of elastomers, there is a continuous 
distribution p(r) of scales r. 

In the gel phase (i.e. r < 0), we find that Eq. I|9a|l solves Eq. Q provided that 



2H 



and that p(r) satisfies the following integro-differential equation: 



(10) 



-p'(r) 



4w 



p(r) 



4 v 



p(t')p(t — t ) dr' . 



(11) 



Equations l|10|) and are identical to those found for the isotropic case of VT [12]. The 
stability of this nematic saddle point can also be established in a way similar to the isotropic 
case (see Ref. [20]). 

We now come to the main point of this Letter: obtaining the elastic free energy of isotropic 
and nematic random solids. To do this, we shall impose an arbitrary homogeneous deformation 
of the boundary of the system, which we encode in the matrix A and illustrate in Fig. ^ We 
shall not make any assumptions about how the interior of the system changes in response to 
this deformation. Our aim is to determine the new saddle point Cl\ that minimises the free 
energy and is consistent with the deformation of the system. (We explain this consistency 
further, below.) We proceed by hypothesizing a modification of the original saddle-point 
solution ij^ajl: 



x — z 
1 



x" - A • z 
1 



(a= 1,.. .,n), 



l+n 



V V" 



(12a) 
(12b) 



In general, det A (= V/Vq) may differ from unity, corresponding to a change in system volume. 
By substituting this modified Ansatz into Eq. l(7)l. we find that the Ansatz is indeed a solution, 
provided p(r) is the same distribution as defined by Eq. (|11|) (in the limit n — > 0). 



( 5 )This automatically ensures that the measurement ensemble has the same volume as the preparation en- 



semble, i.e. V 



Vo. 
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The interpretation of this saddle point is as follows. After the deformation, the same frac- 
tion of the monomers (i.e. q per unit volume) are localized. In the preparation ensemble, a 
localized monomer continues to exhibit Gaussian fluctuations around the point z with an un- 
changed variance-matrix t _1 1o. However, in the measurement ensemble it fluctuates around 
the new point A ■ z but with the original variance-matrix r _1 l. This implies that the average 
position of each monomer, parametrized by z, is deformed afjinely, whereas the fluctuations 
around z remain intact ( 6 ). Note that z, as well as x , are confined to the volume of the 
preparation ensemble (i.e. the range over which z is integrated; see Fig.0. Observe, further- 
more, that Cl^ vanishes whenever x" and A • x° are widely separated (for any a = 1, . . . , n). It 
follows that the x" are confined to the preparation ensemble transformed by the distortion A 
(i.e. the measurement ensemble), and this establishes that A is the homogeneous deformation 
imposed on the boundary of the system (see Fig. QJ. ( 7 ) As we shall show in a forthcoming 
paper [22], the affine character of the deformation is destroyed by the fluctuations of nematic 
order Q° in the preparation ensemble that become frozen in at the time of cross-linking. 

We now calculate the elastic free energy density of nematic elastomers / e j(A) at the 
mean- field level. To do this, we insert the deformed saddle point l|9a|l . modified according 
to Eqs. (|12|l . into the total free energy density (JSJ), and subtract its value for the undeformed 
(A = I)saddle point. Then, dividing appropriately by n (recall that there are n replicas of 
the measurement ensemble), and taking the replica limit n — > 0, we find 

/ei (A) = Umi(ff[n A ]-£r[fi]) 

= A iTrl A T r 1 A + — — -B(detA-l) 2 , (13a) 

» = ~\rf, B = B 0P 2 . (13b) 

It is clear that the first term in Eq. I|13a|) coincides with the free energy density of the neo- 
classical theory of nematic elastomers, Eq. (J2J; the second term describes the energy cost for 
volume changes, with the bulk modulus B related to the parameter Bo, Eq. (|13b() . Con- 
sequently, what we have derived, Eqs. I|13|l . is the neo-classical elasticity model of nematic 
elastomers, in fact generalized to finite bulk moduli systems. In the limit B — > oo, the incom- 
pressibility constraint det A = 1 is restored. Furthermore, if Q° = Q = we have lo oc 1 oc I, 
and our result trivially reduces to the classical theory of rubber elasticity, Eq. QJ. Finally, 
that the shear modulus, Eq. l|13b|) . scales as |r| 3 is a mean-field result and has also been 
derived via other methods [21,23]. 

We emphasize that Eq. (|13afl is derived from the Landau theory of VT, which includes 
the the most relevant contributions. Therefore it is independent of short-distance details, and 
thus provides a universal mean-field description for the elasticity of all forms of vulcanized 
matter near the vulcanization point, provided that the corresponding transition is described 
by the Landau theory. This observation may explain, in part, the huge success of the classical 
theory of rubber elasticity, Eq. JIJ, an d its anisotropic generalization, the neo-classical theory, 
Eq. ©. 

The present work also constitutes a starting point for studying spatial fluctuations, both 
thermal and quenched, in vulcanized matter of various forms [24]. Sufficiently close to the 

( 6 ) A similar result was found for isotropic VT systems [21]. 

( 7 ) For careful readers we note that, in the isotropic limit, this deformed saddle-point solution is slightly 
different from the types of Goldstone fluctuations studied in Rcf. [23], partly due to different parameteriza- 
tions. A detailed discussion of the connections and differences between these two analyses will be presented 
elsewhere [19]. 
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vulcanization point, critical fluctuations of the VT order parameter field f2 become important, 
and change qualitatively the scaling of fj,. (This issue was recently addressed; see Ref. [16].) 
Nevertheless, the form of the elasticity theory J5J continues to hold, upon the incorporation 
of critical fluctuations. On the other hand, fluctuations of Q° in the preparation ensemble — 
which may be strong at short length-scales even deep in the isotropic phase — would provide 
a source for quenched random stresses in the random solid state, and these would couple 
directly to strain. Consequently, quenched-in fluctuations in Q°, and their alignment by 
external stress, may change the elasticity qualitatively, and are likely to be responsible for 
the universal deviation of the stress-strain curves from the classical theory, observed for half 
a century. This issue will be explored in forthcoming work [22] . 
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